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Abstract 



We consider the relation of the multi-component 2D Toda hierarchy with matrix orthog- 
onal and biorthogonal polynomials. The multi-graded Hankel reduction of this hierarchy is 
considered and the corresponding generalized matrix orthogonal polynomials arc studied. In 
particular for these polynomials we consider the recursion relations, and for rank one weights 
its relation with multiple orthogonal polynomials of mixed type with a type II normalization 
and the corresponding link with a Riemann-Hilbert problem. 

1 Introduction 



' Recently, multiple orthogonal polynomials, the related Riemann-Hilbert problems and its ap- 

plication to different areas, for example Brownian motions have deserved much attention pQ. 
The field of matrix orthogonal polynomials has been also a growing area of research with some 
similarities with the scalar case but much more richness, see [2]. The relation of multiple orthog- 
onal polynomials with multicomponent KP hierarchies has been noticed in [3] and the string 
equation formalism of integrable systems has been applied in Ten years ago Adler and 
van Moerbeke [5] -in the context of the so called discrete KP hierarchy- introduced what they 
named as generalized orthogonal polynomials, and what they claimed to be the corresponding 
Riemann-Hilbert problem, while later on they studied the related Darboux transformations [6] . 
Recently, and for the Toeplitz case Cafasso [7j extended this work in order to consider block 
matrices within the non-Abelian Ablotiwz-Ladick latttice. 

Following Ueno and Takasaki [8] and the seminal paper of Mulase [9] in [10] we gave a 
description of the infinite multicomponent 2D Toda lattice hierarchy in terms of a Gaussian 
factorization (also known as Borel factorization) in an infinite dimensional Lie group, while later 
on [11] we analyzed the dispersionless limit of the hierarchy with the aid of the factorization 
problem. (See [12] for a discussion on different cases where this factorization makes sense.) 
Following [5] we could argue as follows: i) in the one hand the multicomponent Toda hierarchy 
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may be viewed as an LU factorization of certain deformed infinite-dimensional matrix and ii) on 
the other hand, the same matrix could be thought as a moment matrix and the corresponding 
LU factorization should give us the corresponding generalized matrix orthogonal polynomials. 
In this manner, we would be able to build a bridge between multicomponent Toda hierarchy 
and matrix orthogonal polynomials. This is the main idea developed in this paper. First, 
we connect matrix orthogonal and biorthogonal polynomials with the multi-component 2DToda 
lattice hierarchy, focusing in particular on the Hankel reduction. Second, we generalize the band 
condition of [5] to the multicomponent case and consider what we refer as multi-graded Hankel. 
This leads to a multi-component extension of generalized orthogonal polynomials which in some 
cases can be described in terms of multiple orthogonal polynomials of mixed type with a type II 
normalization. This connection allows us to give an appropriate Riemann-Hilbert problem for 
this generalized orthogonal polynomials (notice that the one discussed in [5] is not correct). 

The layout of the paper is as follows. In this introduction we give a overview of the Gaussian 
factorization and the semi-infinite multi-component 2D Toda lattice hierarchy. Then, in §2 we 
discuss how matrix orthogonal polynomials and the multi-component 2D Toda lattice hierarchy 
are connected. Finally, in §3 we consider multi-graded Hankel reduction, extended generalized 
orthogonal polynomials, mixed multiple orthogonal polynomials and corresponding Riemann- 
Hilbert problems. 

1.1 Gaussian factorization and the semi- infinite multi-component 2D Toda 
lattice hierarchy 

For the construction of a Lie group theoretical setting we denote by A the shift operator for 
matrix valued sequences. The associative algebra of linear operators on these sequences can 
be identified with the associative algebra of semi-infinite matrices with entries taking values in 
C NxN , the set of N x N complex matrices. With the usual commutator for linear operators this 
algebra is also a Lie algebra denoted by whose Lie group G is the group of invertible linear 
operators in q. Let's take g £ G and consider the following Gaussian factorization problem 
g = S~ l S where S is a block lower triangular matrix, with Su = In, being I/y £ £, NxN the 
identity matrix, and S is an block upper triangular matrix. In [12] it is proven that the Borel 
decomposition holds if all the principal minors do not vanish. Thus, the factorization holds 
under "small" continuous deformations and we can consider the factorization g(t) = S(t) S(t) 
where t stands for a set of complex variables. As was discussed in [TO] this factorization problem 
leads to an integrable hierarchy of nonlinear PDE known as multicomponent 2D Toda lattice 
hierarchy. Let us discuss these issues in more depth. Observe that the matrix associated with the 
shift operator is the block matrix (IjvA,i+i) an d A* is the operator associated with the transposed 
matrix (Ijv^i+i,i)- If E a b, a,b = 1, N is the canonical basis of C NxN and t = ({tj a }, {ij a }), 
j = 1,2,... and a = 1, . . . , N, is a collection of complex parameters we introduce 



A' 



W 



J2E aa expC^2t ja A j ) 



0=1 j=l 



N 



W : 



^K a exp(^t- a (A*) J ) 



o=l 3=1 



and consider Gaussian factorization of g(t) := Wo(t)gWo(t) 



-i 
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Following [TO] we define the Lax operators 



L 


:= SAS- 1 = 


A + u + u l A t + u 2 {A t ) 2 + 


C a 


:= SE aa S 


= E aa + c a iA* + ■ • • , 


L 


:= SA^- 1 = 


= e*A i + u + niA + -- - , 


C a 


'■= SE aa S 1 


= C a0 + C a lA H , 



(1) 



where all the coefficients in the A-expansions belong to C NxN . The multi-component 2D Toda 
hierarchy has the following Lax representation 

<>L [&C a )+,L], ^ = [(UC a ) + ,C b ], ^L = [(VC a ) + ,L], ^ = [(lJC a ) + ,C b ], 



ULja UL ja ub ja uo ja 

c)T BC f)T Be 1 

^ [{UOaULl b={(LiC a )-,C b }, — = [(Z^7 a )_,Z], <L = [(lJC a UC b ], 



ULja UL ja ub ja ub ja 

where the sub-indices + and — denote the block upper triangular, strictly block lower triangular 
projections, respectively. 

2 Matrix orthogonal polynomials and the multi-component 2DToda 
lattice hierarchy 

Following Adler and van Moerbeke [5] we construct families of matrix orthogonal and bi- 
orthogonal polynomials associated with the 2D Toda lattice hierarchy. 

In the first place, we define the following families of (time-dependent) matrix polynomials 

PO) = {Pi(z)}i^ := Sx{z), p{z) = {pi(z)}i^o ■= (S^xiz), 

where x( z ) := Obv> zl^, z 2 In , . . .)* and the symbol ' denotes Hermitian conjugation. Next we 
consider a matrix-valued bilinear pairing between matrix polynomials. Given matrix polynomials 
P{ z ) = Ylk=o PkZ k and Q(z) = ^2 J i =0 Qiz l (of degrees i,j, respectively) we have 

(P(z),Q(z)) = Y, Pk(z k I N ,z l I N )Ql 

k=l,...,i 
1=1,.. .J 

where (z k lN, z 1 In) denotes the matrix for the bilinear pairing in the canonical basis and for each 
(k, I) is an N x N complex matrix. 

This pairing has the following properties 

1. Is linear in the first component: 

(ciPi(z) + c 2 P 2 (z), Q(z)) = Cl (Pi(z), Q(z)) + c 2 (P 2 (z), Q(z)), Vci, c 2 G C NxN 

2. Is skew-linear in the second component: 

(P(z),c 1 Q l (z)+c 2 Q 2 (z)) = (P(z),Q 1 (z))c\ + (P(z),Q 2 (z))cl Vc x , c 2 G C NxN 
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Proposition 1. 1. If (2'I/v, z^n) = 9(t)ij where gij is the C NxN block in the position (i,j), 
then the families p(z) and p(z) are biorthogonal matrix polynomials for the linear pairing, 
i.e. 

{Pi(z),Pj(z)} = 5ijI N . 

Moreover, 

( Pl (z),z l I N ) = 0, l = 0,...,i-l, (z l l N ,Pj{z)) = 0, Z = 0,...,j-1. (2) 

2. In addition, if the time- dependent initial condition g(t) is Hermitian for all t then p(z) 
and p(z) are two families of matrix orthogonal polynomials, moreover, the two families are 
proportional. 

Proof. 1. With the previous definitions for p{z) and p(z) we have: 

P i(z) = j2s ik z k , Pj (z) = j2(sr^z l , 

k=0 1=0 

where and S^ 1 are the blocks (i, k) and for S and <S _1 , respectively. Hence 
(pt(z),Pj^)) = E S ik (z k l N ,z l l N )S^ = Y, S ik {z k l N ,z l l N )S^ 

k,l=0 k,l^0 

= (Sg(t)S-% = (SS-'SS- 1 )^ = SijlN, 

as desired. Finally, ([2]) is proven by induction. First we have that (pi(z),po(z)} = 
0, but po(z) = (Sqq)^ is invertible and therefore we conclude Ijv) = 0. Now, 

(pi(z),pi(z)) = 0, but pi(z) = (S^y z + (S^q)^ , and using the skew-linearity, the previous 
result and the fact that (^n 1 )^ is invertible we deduce that (pi(z), zIn) = 0, and so forth 
an so on. 

2. Let's study the conditions under which p(z) is a family of matrix orthogonal polynomials. If 
we take two polynomials in the family, such as Pi{z) = J2k=i SikZ k and Pj(z) = J2j=i Sjiz 1 
we have 

i,j oo 

{Pi(z),Pj(z)) = Y S lk (z k I,z l I){S^ = Y, S ik (z k I N ,z l I N )(S% 

k,l=l k,l=0 

= {Sg®^ = (SS%. 

Observe that SS^ is clearly block upper-diagonal with its Hermitian conjugate given by 

(S5 t ) t = S& = Sgftigit)- 1 )^ = S^i)" 1 ^ = SSl 

Therefore, SS' is Hermitian and block upper-diagonal, which implies that SS^ is a block 
diagonal matrix and the blocks in the diagonal are <C NxN Hermitian matrices. 

We conclude that (pi(z),pj(z)) = 5ij{hi)~ l , where hi is a Hermitian matrix. Notice also 
that as a consequence p{z) = hp(z) where h = diag(/ii, /12, . . .). 

□ 
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2.1 The Hankel case 

We choose g to be a block Hankel matrix so that Ag = gA l or 

(Ag)ij = g i+ld = g iJ+1 = (#A%-, 
In this case we have for the blocks of the moment matrix g 

for some matrices 7^') G C NxN . From 

Ag(t) = AWogWu 1 = W^AgW^ 1 

= WogA'W^ 1 = WogW^A* = g(t)A l , 

we easily deduce that g(t) is block Hankel if g is. We have 

Proposition 2. Assume that g is block Hankel then with j^) a block moment matrix; i.e. 
jti) = J R xi p(x)dx then the pairing can be viewed as a scalar product in the real line whose 
matrix moment is g ; i.e., 

(P(x),Q(x)) = [ P(x)p(x)Q(x^dx. 

Proof. On one hand we have (P(x),Q(x)) = Pil^^Q]- Using the previous definition 
f R x : > +k p(x)dx = (x^Im , x h lN) = j(j+ k ) (the Hankel symmetry ensures that there is only depen- 
dence in j + k) we have: 

(P{xlQ(x)) = Y,Pil {i+j) Q) = Y. P i J x 3+k p(x)dxQ} 

= f P(x)p(x)Q(x^dx. 
Jr 

□ 

In general arbitrary continuous deformations do not preserve the Hermitian character of g. If 
we look for families of matrix orthogonal polynomials on the real line we should make restricted 
deformations. Let's make this point clear. In the following z* denotes the complex conjugate of 
z G C. 

Proposition 3. If the matrix g is block Hankel and the matrices 7^ are Hermitian then 

1. The families p(z) and p(z) are proportional and the two of them are matrix orthogonal 
polynomials in the real line. 

2. Moreover, if the continuous deformation parameters satisfy one of the two following con- 
ditions 

(a) tj a , tj a G R and satisfy tj a = tj,ij a = ij, a = 1, . . . , N. 

(b) t ja , ij a satisfy t ja + i* a = 0. 
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the result holds for the time dependent moment matrix. 

Proof. 1. If the matrix g is block Hankel and the blocks are Hermitian the matrix g is itself 
Hermitian. Given pair of indices for an element of g, there exist four integer indices 
(k,l),(m,n) with k, I ^ y m,n = 1, . . . ,N that satisfy = (A M ) mn = (A k i)* nm = 
(Aik)nm = a ji ( we use A € C nxN for a block of g) as a consequence g = g^ and by 
Proposition [JJ the first part of the result holds. 

2. Let us compute the time evolution for g. We have 

N oo N oo 

Wo = ^ s «» ex P£ *i« Ai ) = E ^ E s f A "> 

a=l j=l a=l j=0 

where is the j-th Schui0 polynomial for the component a. Consequently and taking 

oo N oo 

j=0 k=l j=0 

so it is straight forward to see the block structure of Wq, whose blocks are given by 
(Wo)ij = s j-i if J — z ^ and Oat otherwise. 

A similar argument for Wq 1 leads to (Wq )jj = Sj_j if i — j ^ and Oat otherwise (here 
Sj := X^aLi Ea, a Sj and s*- is the j-th Schur polynomial for the component a but now in 
the variables —ij a ). 

We are now ready to compute 

<?(% = £ (Wok^^o -1 )*; = E *k-*7 (fc+0 sj-j = E »fc7 (<+i+fc+,) 5i. 

fc,Z>0 k^id^j k,l^Q 

Then 

(a) If the first condition holds then all matrices Sj and Sj are real and scalar so 

A:,/>0 fc,Z^0 fc,Z^0 

= g(t)ij. 

(b) If the second condition holds sj = sj for all j so 

g (t% = J2 sf 7 (i+i+fc+,) 4 = E 47 (i+j ' + ' +fe) 4 = E «*7 (<+i+fc+,) si 

fc,/>0 Z,fc^0 fc.^O 



1 We remind the reader that the Schur polynomials are given by exp(^°l 1 tj a h?) = X^jlo s j ^ so s o ~ 1 
and sj a ' = J3p=i Sji-i hj =» * J i a ' ' ' * J p a ~ * J£I + ' ' ' wriere the missing term is a polynomial of utmost degree j 
in the variables tio, *2a, • - ■ , tj-i,a- 



G 



Under any of the two conditions g(t) is block-Hankel with Hermitian blocks, so is itself 
Hermitian. That proves the second part of the proposition. 

□ 

We now discuss the recursion formulae. Using the Hankel condition Ag = gA l and taking the 
usual definition for the Lax operator we conclude L = L and hence we have L = A+u(n)+i> (n)A t . 

Proposition 4. The polynomials p(z),p(z) satisfy a three term recurrence law given byp n+ \(z) = 
zp n (z) - u(n)p n (z) - v(n)p n -i(z) 

Proof. Prom the definition of L we have Lp(z) = LSx(z) = zp(z) = (A + u(n) + v(n)A t )p(z). 
If we take the sequence terms {p n (z)} n we conclude zp n (z) = p n+ \{z) + u(n)p n (z) + v(n)p n -i. 
For the polynomials p(z) we have L^p(z) = zp{z) = (Aim (n) + u^(n) + A t )p(z), and using the 
sequence {p n {z)} n we obtain another recurrence law. □ 

3 Multi-graded Hankel reduction, generalized orthogonality, mul- 
tiple orthogonal polynomials and Riemann— Hilbert problems 

In this section we will study generalized Hankel type conditions. Given a multi-index n = 
(m, . . . ,njv) with n a non-negative integers we define for A G g the power A n = Yla=i A na E aa . 
For two multi-indices n and m a matrix g is said to be a (ft, m) multi-graded Hankel if 

A n g = g(A t )" 1 . (3) 

If as before gij G C NxN denotes a block in g then we can write gij = (gij : ab)i^a,b^N and the 
multi-graded Hankel condition reads gi+ na j,ab = 9ij+m b ,ab- An ample family of multi-graded 
Hankel matrices can be constructed in terms of weights Pj, a b as the moments 



9ij,ab = / x l p jt ab(x)dx, (4) 
where the weights satisfy a generalized periodicity condition of the form 

Pj+m b ,ab(x) = X na p jtab (x). (5) 

Thus, given the weights po, a b, ■ ■ ■ ,Pm b -i,ab, all the others are fixed by ([5]). From now on, we 
concentrate only in these cases of multigraded Hankel matrices. We notice that for the 1- 
component case and for n\ = mi case these moments matrices were studied in [5]-[6]. 

Proposition 5. For multi-graded Hankel matrices the matrix polynomials pi satisfy the following 
generalized orthogonality conditions 



J 

Jr 



p i (x)p j (x)dx = 0, j = 0,...,i-l, Pr-=(pj,ab)£C NxN . (6) 

Proof. From Sg^S^ 1 = 1 G we get ^ j=i,...,i Pij,ab9jl,bc = for a, c = 1, . . . , N and l = 0,...,i- 

b=l,...,N 

1. Now recalling ^) we get the result. □ 
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Using the Euclidean division i = 6 c m c + a c , with 9 C ^ 0, ^ a c < m c we get a better insight 
of the orthogonality relations ©, for a, c = 1, . . . , N, 



N 

V] / p l ,ab(x)p j .bc(x)(x nb ) l dx = 0, j = 0,...,m c -l, Z = O,...,0 C -1, 



6=1 

iV 



/ Pi,afe(x) / o i , bc (x)(x ni ') 9c dx = 0, j = 0, . . . , a c - 1. 



(7) 



3.1 Evolution 

From @ we conclude that is of (n, m) multi-graded Hankel type if g is. In general, the 
evolution of g is given in terms of Schur polynomials by 

9ij(t) = S k9i+k,j+lSl 

which recalling leads to the following evolution of the weights 

P&) = ^2 xks kPj+isi = exp(t(x))^2 p j+ isi, 



(8) 



k,l>0 



l>0 



where t(x) = J2 a =l ^a{x)E aa and t a (x) := tjaX 3 ■ It can be easily checked that Pj(t) satisfies 

the periodicity condition © if pj does. From ([8]) we infer that 



/ f>0,ab(t) \ 

: = exp(t a (x)) 


( ^0,ab 
X " ' ■S^ra b — \^ab 
X °^m|,-2,a() 


^l,ab 
^0,ab 
X a ■5^m b ~\,ab 


^ 2,o6 
^ D,o6 


'5^'m b — l,ab^ 

^m b —2,ab 

^m b —2,ab 


\Pm b ~l,ab{t) J 


V x n ^ 1:ab 


X n ^2,ab 






where 












y= 1 


mj— 1 

£ 4 fe exp(- 

fc=0 


-t 6 (£fc"° /m6 )) 


E b — 




If we denote 














/ = 0, 1, . . . ,m b - 1, 


we have 










t a (e k b x) = i^(x)+^(x), 


f ] (*) = 


=l,...,m&~l 


kl jm b +l 
-l £ b x ' 



/ PO. 



ab 



\Pm b — l,ab / 

(9) 



1. 
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and therefore 

m b -l 

^ = ^W^ eX P(-4 0) (^ /mt )) £ 4 fc exp(-tf (x-/-)). (10) 

k=0 

Finally, we deduce 

m b -l 

Phab (t) = eMta(x) - (x^ m »)) /5^exp(-f ] (x^)), (11) 

k=0 

where we have used the discrete Fourier transform of the weights 

m b -l 

3.2 Recursion relations and symmetries 

In terms of Lax operators the multi-graded Hankel reduction reads |10j 

N N 

£>.-J2L n «C a = ^2L mb C b . (12) 

a=l 6=1 

Within this subsection we assume that 

m ^ • • • ^ njv ^ 1, 77ii ^ • • • ^ m-Tv ^ 1, 

and suppose that n\ = ■ ■ ■ = n r and n r > n r+ i. Given (JTJ) from (|12p we deduce that 

JSf = (En + • • • + E rr )A ni + jSf ftl _ 1 A ni - 1 + • • • + JSf + J2LiA* + i?- 2 (A*) 2 + • • • , 
while we also have «Sf = X^tLi L mb Cb and therefore 

JS? = JSf_ mi (A*) mi + JSf-^+iCA*)" 11 " 1 + • • • + J% + JSfiA + if 2 A 2 + • • • . 
We conclude the block band structure 

if = (E X1 + ■■■ + E rr )A n ^ + J^-iA" 1 " 1 + • • • + jSf^+xCA*)™ 1 " 1 + if_ mi (A*) mi . (13) 
Proposition 6. T/ie polynomials Pi(z) are subject to 

n 

(E u + ■■■ + E rr ) Pi+ni (z) + • • • + if_ mi pi_ mi (z) = pi(z) ( ^ z na E aa ) 

a=l 

Proof. We only have to show that J£p = p(^2™=i zUa E a a)- But this follows from Jfp = 

Similarly, from ££^p{z) = p{z)(J2bLi zmb ^bb) a recursion relation follows for the p{. 
Finally we notice that in this case the following symmetry conditions hold [10] 

N d N d N d N d 

a=l ,n «° a =l aa x a=l ul m a a a _ l < J^im a a / 

N a N a N a N a 

E 5 i 2 -+E5r-)Q = (E5r-+Ear-)ft = * 

for i 1 and 6 = 1, . . . , N. 
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3.3 Relation with multiple orthogonal polynomials 

We will see that when the weights pj are particular rank one matrices there is a nice corre- 
spondence with multiple orthogonal polynomial of mixed type with a normalization of type II. 
Following [1\ we take two sets of non negative multi-indices v = (ui,..., v p \ = (tii, . . . , fj, q ) 
and write \u\ = Yl P j=i U J an d H = Y^k=\1 1 k- We also take weights {wu} p J=1 and {w2k}k=i 
which are assumed to be non-negative functions on the real line. For a fixed pair u, ft we say that 
{^p,p,,j}j=i,...,p is a set of multiple orthogonal polynomials of mixed type if degAp^j ^ vj — 1 
and the following orthogonality relations 

y^Ap tfi j(x)w 1 j(x)w2K(x)x a dx = 0, a = 0, ... , /jl k ~ 1, K = l,...,q, (14) 

are satisfied. Alternatively defining the following linear forms Qp,p,{x) := ]Cj=i J ^u,p,,j{x)wu(x) 
the orthogonality relations can be written in the following way 

/ Qp,fi(x)w2K(x)x a dx = 0, a = 0, ...,/XK-l, K = l,...,q. 

Jr 

For multiple orthogonal polynomials as described above we will take \v\ = and assume that 
the following conditions hold 

1. For each K = 1, . . . ,p the orthogonality relations for the multi- indices v + e^,/i have a 
unique solution with A d ^^k monic and degA^^^x = — 1, with deg^lp^j < vj — 1 
if J 7^ if. We will call it type II normalization to the if-th component and write that 
normalized solution as {^4?^ j }j=i,...,p- 

2. For each if = 1, ... ,q the orthogonality relations for the multi-indices v,p, — ex have a 
unique solution with the following normalization: deg ^4p,^,j = v j—1 and L Q u ^{x)w2KX^ K dx 
1. We will call it normalization of type I to the K-th component and write that normalized 
solution as {A { *£)}j= lt ... >p . 

In order to connect ([7]) with multiple orthogonal polynomials we consider the Euclidean division 
i = Qan a + ?"a,with q a ^ and ^ r a < n a , and write 

n b 

Pi,ab( z ) =J2 zj ~ 1U vMz nb )> (15) 



where IJy- a i,(z nil ) are polynomials in 2™ 6 such that 
deg U ij:ab ^ 



j < r& or j = r& + 1 and b = a, 
1; j > r fe + 1 or j = r;, + 1 and b ^ a. 



Notice that the monic character of pi gives the normalization of Ilj ra+ i aa which happens to be 
a monic polynomial with degllj ra +i,aa = qa- 

The inversion formula for ([To]) can be deduced as follows. If we denote by e b := exp(2-7ri/nb), 
a primitive n^-th root of the unity, and evaluate at e£z, k = 0, . . . , n b — 1 we get the following 
system of equations 

n b 

Pi,ab( e b z ) = J2^b z y~ 1U idM zn ") k = 0,...,n b -l, 
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that we solve in order to obtain the polynomials Hik,ab in terms of a discrete Fourier transform 
of the polynomial through the formula 



fe=0 



Then, ([7|) can be written as 

N n b 



/ ^^n^ ab (x n6 )x J -VMc^)(^ ni ')'dx = 0, fc = 0,...,m c -l, Z = 0,...A-1, 

6=1 i=1 (16) 

„ N n b y > 

/ EE n ^( xnt ) xi "V^)(^ n,, ) ec dx = 0, k = 0,...,a c -l. 
J R b=i j=i 

These equations strongly suggest to perform a change of variables in each integrand of the type 
y _ x n b _ p Qr ^h a t aim, it is relevant that when n& is an even number then supp(/9 j-& c ) C M + , 
otherwise the change of variables is ill defined. In fact, for these even cases, one easily see that 
the weights must be supported on the positive axis or uniqueness of the orthogonal polynomials 
is not ensured. Moreover, for n\, odd it is also necessary to assume that the weight is supported 
either only in the positive real numbers or only in the negative real line, this requirement comes 
from the positivity condition on the weights and the use of (|19|) . Hereon we will assume that the 
all weights are supported on the positive real semiline. When the mentioned change of variable 
is performed in (|16p we get 



N n b 



/EE u ijMy)Pjk,bc(y)y ld y = °i fc = o, . . . , m c - i, i = o, . . . , e c - i, 

J R 6=1 3=1 
,. N n b 

/EE u ij,ab(y)pjk,bc(y)y dcd y = °' k = o,...,a c -i, 

Jr 6=1 7 =1 



(17) 



with 



PjkMy) = —y" b Pk,bc{y nb )■ 

n b 

Now, if the matrix weights pk are rank-1 matrices of the following particular form 

Pk,bc( x ) = vi, b {x)w 2) kc{x nb ) 1 

we get 



,. N n b 

/EE u ij,ab(y) w ijb(y)w 2 ,kc(y)y ld y = °> = o, . . . , m c - i, i = o,...,e c -i, 

J R 6=1 3=1 

„ N n b 

/EE n ^>6(y)^i ,jb(y)w2,kc(y)y 6cd y = °; fc = 0, . . . , <T C - 1, 

^ M 6=1 3=1 



(18) 
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where 



1 J--i JL 

wijbiv) = — V nb v hb (y"b). (19) 
n b 

We are ready to describe the relation among multiple orthogonal polynomials of type II and 
generalized matrix polynomials. First, given (a, j) with a = 1, . . . ,N and j € N we make the 
definitions ftf(a,j) := n\ + ■ ■ ■ + n a _i + j and M.(a, j) := mi + • • • + n a -i + j + 1. 

Proposition 7. Relations (|18j) are particular cases of (114p mi/i: 

1. J = Af(b,j) for b = 1, . . . ,N and j = 1, . . . , ra;,; and = .M(c, /c) /or c = 1, . . . , N and 
k = 0, . . . , m c — 1. I^e Ziawe therefore the identification J = n\ + • • • + n^-i + j and 
K = mi + • • • + m c _i + k + 1. 

£ p = |n| = m H h ?2 7v and q = |m| = mi + ■ ■ ■ + m^. 

J ^ r b or j = r b + 1 and b = a, 
J > r b + 1 or j = r b + 1 and b ^ a, 

o~ c ^ k ^ m c — 1. 



3. 



^M{c,k) 



Qb + 1, 
ft, 



^- 1^1 = = -Wi + 1 <md = YIk = Ni. 

5- A v ^j = U ijtab andp i>ab (z) = zJ ~ lA u,fi,M(b,j)( znb )- 

6. degAp^j = uj if J = Af(a,r a + 1) so that Ay^j = j^ I _^ < - a ' ra+1 ^ ; an d we have a type II 
normalization w.r.t J\f(a,r a + 1). Thus, Pi A b{z) = Y^L\Z j ~ 1 A^^^ +1 '\z 1lb ). 

The reader should notice that evolution of the weights given trough ([8]) give the following 
evolution of wij b and W2,k 



wi,jb{t) = exp(E y k/nb t kb )wij b , W 2 ,kc{t) = E W 2,k+l cs\ 



k>l l>0 



Using (llip we may write 
w 1>jb (t) = ex.p(^2y k/nb t kb )w 1Jb , 



m b -l m b -l 

,r 2M (t) = eM-^ (y 1/mc )) E 4tex P (-tf(y 1/mc )), <L : = — E ^y~ t/mc ^k +i , c 

1=0 



These evolved weights fulfill the positivity condition (recall that their support is included in the 
positive real line) when we have tj b ,ij mc c € K and ij c ^ for j ^ mod(m c ). 
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3.4 Riemann— Hilbert problems 



As we have just shown the generalized matrix polynomials are connected with a family of multiple 
orthogonal polynomials for a particular rank one moment matrix. In [1] the Riemann-Hilbert 
problem for multiple orthogonal polynomials of mixed type was presented, see also |13| . We will 
discuss its relation with the generalized orthogonal polynomials Pi(z). 
Let us recall the reader that the Cauchy transform is defined by 



Qv',p,',K{z) := -— / - W2K{X)dx. 

* m JR Z — X 



Now we can make the following definition for the (p + q) x (p + q) complex valued matrix 
Y(z) 

Yiri — A (II ' K) 
v ._ A( n > K ) 

*K,J+p ■— W D ' +eK ^j 

Yk+ p ,j ■= -2TTiA { J f 'p_ eKtJ 

Y K +p,J+p ■= -^Q { ^- eK j 



J = 


1,.. 


• ,p 


K = 


1,. 


• • ,P 


J = 


1,- 


■,Q 


K = 


1,. 


■ ■ ,P 


J = 


1,- 


• ,P 


K = 


1,- 


..,q 


J = 


1,. 


■ ,q 


K = 


1,- 


..,q 



we will also use the following real valued (p + q) x (p + q) matrix D(x) defined by blocks 



D(x) 



Hp W(x) 



qxp 



where Wjk(x) = wij(x)w 2 k(x). 

We adapt to the present situation a result of pQ taking into account the support of the 
weights 

Theorem 1. Let be v',fi' two multi-indices such that \v'\ = \p!\and suppose that the normality 
conditions hold. Let also be two sets of weight functions {wu}j=i t ,,, :P and {u)2K}K=i,...,q such 
that for every J,K wu, (x)w2k(%) &re differ entiable a.e. in M+ and x^wu,x^W2K £ 
j = 0, . . . , v' K — 1. In x = we require the weight functions to be bounded. Then, the matrix 
Y(z) is the only solution of the following Riemann-Hilbert problem. 

1. Y(z) is analytic in C \M+. 

2. Y(x) + = Y(x)-D(x) for all x > 0. 

3. Y(z) diag(z-^, . . . , z-4,*"v . . , z^v) = I (p+q) + Q for z 



oo. 



I Y(z) = O 



( 1 1 
1 1 



1 log \z\ 
1 lofflzl 



log I z I 
log I z I 



1 log|z| log \z\ 



log \z\\ 
log \z 

log |z|/ 



when z — > 0. 
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Proof. First we see uniqueness of the solutions. First we see that as det-D(x) = 1 det Y(z) 
is analytical across the integration contour, then the only possible singularity of det Y(z) is in 
z = 0. As zdetY(z) — > when z — > the isolated singularity must be removable and det Y(z) 
is an entire function. The asymptotic behavior and Liouville's Theorem gives detY(z) = 1. 
Consequently Y(z)^ 1 exists. Given two possible solutions for the problem Y(z) and Y(z), the 
matrix Y(z)Y(z)~ 1 can be singular only in z = 0. As before the singularity must be removable 
so Y(z)Y(z)~ 1 is entire and bounded, and hence Y(z)Y(z)' 1 = I p + q - 

Now we prove that the matrix Y(z) is a solution of the R— H problem. Condition Q] fol- 
lows from the fact that polynomials are entire functions and general theory of Cauchy integrals 
[14) gives analytic behavior outside the integration contour. Now, observe that condition [2] 
is a jumping condition on the positive real axis and is a consequence of Plemelj formulae. 
For condition [3] we notice that Y{z)a is a monic polynomial with degree u[ for i = 1, . . . ,p. 
We also see that the leading term of each Y(z)u is for i = p, . . . , (p + q) due to the 

orthogonality relations and the type I normalization. Consequently the diagonal elements of 
Y(z)diag(z- l/ i,...,z~ u p , z^i , . . . , z^p ) are equal to 1 and the rest of them vanish like - asymp- 
totically. Finally, condition 2] is a consequence of the behavior of the polynomials and the Cauchy 
transforms. The boundness condition of the weights at z = makes the Cauchy integrals to 
have at much log-type singularities at z = 0. □ 

The Theorem applies to our situation giving 

Proposition 8. If v' = v — ex{a,r a +l) an d fi' = ft; with v and p, as in Proposition^ we have 

H-ij,ab = Xv(o,r a +l),AA(6j) an d 

Pi,ab{z) 
Y N(a,r a +l),Af(b,j)( znb ) 

We observe that in [5] and [6] these generalized polynomials are considered for the one com- 
ponent case N = 1 and for n = m and therefore fit as a particular example of our polynomials 
and consequently are related with multiple orthogonal polynomials of mixed type. Apparently 
they might noticed this fact later, as in [3] they claim that they considered multiple orthogonal 
polynomials in [5]. However, we must say that the Riemann-Hilbert problem derived there is 
different from the Daems-Kuijlaars problem considered here for multiple orthogonal polynomi- 
als. In fact, the matrix Y considered in [5] is not analytic in the upper-half plane and hence 
fails to satisfy the Riemann-Hilbert problem posed there. 
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= ^ zJ ~ 1 Xtf(a,r a +l),Ar(b,k)(z n ''), 
i=l 

= pi e^^W&K 2 )- 

1=0 



14 



References 

[1] E. Daems and A. B. J. Kuijlaars, J. of Approx. Theory 146 (2007) 91. 

[2] A. Duran and F. A. Griinbaum, J. Compt. Appl. Math. 178 (2005) 169. 

[3] M. Adler, P. van Moerbeke and P. Vanhaecke, Commun. Math. Phys. 286 (2008) 1. 

[4] L. Martinez Alonso and E. Medina, J. Phys. A: Math. Theor. 42 (2009) 205204. 

[5] M. Adler and P. van Moerbeke, Commun. Math. Phys. 207 (1999) 589. 

[6] M. Adler and P. van Moerbeke, Int. Math. Res. Not. 18 (2001) 935. 

[7] M. Cafasso, J. Phys. A: Math. Theor. 42 (2009) 365211. 

[8] K. Ueno and K. Takasaki, Adv. Stud. Pure Math. 4 (1984) 1. 

[9] M. Mulase, Adv. Math. 54 (1984) 57. 
[10] M. Manas, L. Martinez Alonso and C. Alvarez-Fernandez, Inv. Prob. 25 (2009) 065007. 
[11] M. Manas and L. Martinez Alonso, Inv. Prob. to be appear (2009). 
[12] R. Felipe and F. Ongay, Lin. Alg. Appl. 338 (2001) 1. 

[13] W. Van Assche, K. T-R Mc Laughlin, A. B. J. Kuijlaars and M. Vanlessen, Adv. Math. 
188 (2004) 337. 

[14] F.D. Gakhov. Boundary Value Problems, Pergamon Press, New York, (1966). Redited by 
Dover (1990). 



15 



